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Abstract 

Low X deep-inelastic ep scattering data, taken in 1994 at the HI detector at HERA, 
are analysed in the Breit frame of reference. The evolution of the peak and width 
of the current hemisphere fragmentation function is presented as a function of Q 
and compared with e^e~ results at equivalent centre of mass energies. Differences 
between the average charged multiplicity and the multiplicity of e^e~ annihilations 
at low energies are analysed. Invariant energy spectra are compared with MLLA 
predictions. Distributions of multiplicity are presented as functions of Bjorken-x 
and Q^, and KNO scaling is discussed. 
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1 Introduction 



In a previous publication |l|] the HI experiment analysed hadron production in ep deep- 
inelastic scattering (DIS) in the Breit frame 0, Q of reference. We showed similarities 
between hadronic distributions in the 'current' hemisphere of this frame of reference and 
those in a given single hemisphere of an e~^e~ annihilation interaction. A similar anal- 
ysis Q was performed by the ZEUS experiment. The first HI study analysed 1993 data 
resulting from e~p collisions. The 1994 run used incident positrons and this analysis has 
more than six times as many events despite tighter cuts on event quality. Since the 
values investigated are still in general below the kinematic region where anything other 
than virtual photon exchange is relevant, the results of the two analyses are directly 
comparable. 

It is the object of this analysis to study the evolution of the fragmentation function 
and charged hadronic multiplicity distributions. One aim is to continue the comparison 
with equivalent e^e~ results as a test of the universality of properties of the partons found 
within a proton and those of the quarks produced in pairs from the vacuum. Another aim 
is to show sensitivity to the running of the QCD coupling constant and to leading order 
(LO) QCD processes from consideration of such hadronic final states. 

After a description of the relevant parts of the apparatus (section 2) and a discussion 
of the treatment of the data (section 3), this paper describes the Breit frame of reference 
(section 4) and extends (section 5) the earlier analysis of the evolution of the fragmentation 
function with Q^, where = —Q^ is the square of the four-momentum transferred from 
the incident lepton. The average charged multiplicity dependence on Q (section 6) is used 
to study the effects of LO QCD processes. For the first time we present invariant energy 
spectra (section 7) from HERA and HI multiplicity distributions (section 8) in the Breit 
frame. 

2 The HI Detector 

The HI detector is described in detail elsewhere @. Here, we give only a short description 
of the components more relevant to this analysis. We employ a right-handed coordinate 
system with the positive (forward) z axis being in the direction of the incoming proton 
beam. Polar angles are measured from this direction. 

Momentum measurements of charged particles are provided, in the central region of 
the apparatus, by two cylindrical and co-axial drift chambers for (r, 0) measurement 
which have significantly fewer dead cells compared with the situation in 1993. These 
detectors are supplemented by two z-chambers. In the forward (proton) direction the 
equivalent measurements are provided by three Radial and three Planar drift chamber 
modules, the Forward Track Detector All these track detectors are inside a uniform 
1.15 T magnetic field. Track segments from all devices are combined to give efficient 
detection and momentum measurement with 5p/p'^ ^1% /GeV for most of the angular 
range used in this analysis, 10° < 9 < 160°. 

In the polar angle range 4° < ^ < 153° the trackers are surrounded by a fine-grained 
liquid argon (LAr) sampling calorimeter ||^ with lead and steel absorber in the electromag- 
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netic and hadronic sections respectively. The calorimeter cells measure hadronic energy 
flow and the energy of the scattered electron for high events. The LAr calorimeter is 
complemented by a backward electromagnetic, lead-scintillator, calorimeter (BEMC 0) 
covering the angular range 151° < 9 < 176°. Behind this there is a double-layer scintilla- 
tor hodoscope which gives efficient background rejection by time of flight measurement. 
The data are derived from two separate samples corresponding to the scattered positron 
being triggered by, and detected in, either the BEMC or the LAr calorimeter. The tran- 
sition is at a of about 100 GeV^ and these samples are referred to in this paper as, 
respectively, low or high data. The trigger is already fully efficient at positron energies 
of 10 GeV. 



3 Data Selection and Corrections 

The full data sample, triggered as above, consists of some 150 000 events taken at an ep 
centre of mass energy, y^, of 300 GeV. Event kinematic variables used in this analysis 
are calculated using only the scattered lepton, which gives both the best resolution in the 
chosen region and clear freedom from bias on the hadronic system studied. Throughout, 
the events are selected to have an identified scattered positron |]10| with an energy E'^ > 
14 GeV. 

To reject beam associated background, we demand that there be no veto from the 
time of flight system and that a vertex is found for the event within 30 cm of the nominal 
vertex position. The low (high) data referred to in the above section are selected to 
be within the limits 12 < < 100 GeV^ (100 < < 8000 GeV^). The dimensionless 
inelasticity variable, y = Q'^/xs, is required to be in the range 0.05 < y < 0.6. These 
conditions ensure that contamination from mis-identified photoproduction events is below 
the 1% level. Together these cuts remove ~ 45% of the original sample. In addition, a cut 
is placed to safely exclude a further ~ 5% of diffractive events which are not well-modelled 
with the DIS Monte Carlo programs used in this analysis and which have no analogue in 
e+e" interactions. We require that a total cluster energy of at least 0.5 GeV should be 
observed in the region 4.4° < 6 < 15°. 

Events are then selected in which a massless quark would be scattered through 10° < 
6 < 150° using four momentum conservation with the scattered positron. This is inside 
the acceptance of the HI track detectors and thus minimises corrections, at the cost of 
removing a further ~ 24% of events. 

There is a source of error arising from QED radiation which comes about because 
of an incorrect, typically overestimated, value of Q^. The error in the boost leads to a 
miscalculation of the direction of the Breit frame axis. In most cases, this leads to an 
apparently empty, or at least depleted, region of phase space where the scattered quark 
fragments are expected. With improved statistics this has become a significant effect. To 
reduce the size of the necessary QED corrections, the new analysis utilises more severe 
cuts on the mass of the total hadronic system, W, of > 4400 GeV^ both using the 
scattered lepton variables, and separately from the hadronic system variables with the 
Jaquet Blondel method ITT . 



Following all of the selection procedures there are 20 810 events in the low Q sample 
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and 1 250 events at high . 

In addition to these event selections there are also cuts made to reject badly measured 
tracks. Any tracks in the central chambers with transverse momentum below 150 MeV 
and those in the forward track detectors with momentum below 500 MeV are removed. 
We also remove tracks failing minimum requirements on the number of hits in a given 
chamber and the quality of the track fit. In order to have only the primary multiplicity, 
cuts are also made to exclude tracks not originating from the interaction vertex. There 
remains a small excess contribution due to the decays of short-lived strange particles. 

To correct for this excess and for acceptance losses, we utilise the DJANG06 [0 
Monte Carlo event generator. This combines a LO perturbative QCD matrix element 
calculation and the colour dipole model of hadronisation with a calculation of QED ra- 
diative effects. Radiative effects remaining after our event selections have been corrected 
by comparing the results of Monte Carlo calculations with and without the inclusion of 
QED radiative effects. These corrections are at the ~10% level. The detector response is 
simulated using a program based on the GEANT package and the simulated events 
are reconstructed and selected using exactly the same analysis chain as with real data. 
The total bin-by-bin corrections made throughout this analysis are generally well within 
±20% and vary smoothly in any given distribution. 

The largest individual sources of systematic error for this analysis originate in possible 
calibration errors of ±1% for the BEMC and ±3% for the LAr electromagnetic calorimeter 
energy scales. These directly affect the accuracy of Lorentz boosts and give rise to an 
uncertainty in the number of tracks of ~5% at low and ~8% at high Q^, irrespective 
of other kinematic selection. The corresponding systematic error in the Monte Carlo 
derived acceptance correction functions has been estimated to be ~ 2% using several 



different generators |T^, |T^. Visual scans of real and simulated data have ascertained 
that the efficiency of the track detectors is modelled to an accuracy of better than 2% by 
our Monte Carlo simulations. 

The full experimental details of how the corrected distributions of this analysis are 



obtained from the data may be found in reference [ITB] . 



4 The Breit Frame of Reference 

The ep Breit frame is aligned with the hadronic centre of mass (HCM) but boosted along a 
common z direction such that the incident virtual photon has zero energy, zero transverse 
momentum and a z component of momentum —Q. As with the laboratory frame of 
reference, we choose the positive z axis to be in the direction of the incoming proton. 
The negative z direction is referred to as the 'current' hemisphere of the interaction. Our 
earlier Breit frame analysis showed that multiplicities in the current region of the Breit 
frame depend on Q and not on Bjorken x (~ Q^/W^^), as opposed to the HCM where 
multiplicities depend on the natural scale, W In the naive quark parton model 

(QPM) the massless incoming quark has energy Qjl and z component of momentum 
+Q/2, carrying an approximate fraction x of the proton's momentum. After scattering it 
still has energy (5/2, with momentum —Q/2. By comparison, we thus take the equivalent 
to the e+e~ centre of mass energy, E* , to be Q. 
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As compared with the HCM the Breit frame current hemisphere is dominated by the 
fragments of the struck quark alone; the 'spectator' proton remnants go entirely into 
the 'target' hemisphere, with much higher momentum. There is excellent acceptance for 
the current hemisphere in the central HI detector. For example, even in the highest 
interval of this analysis there is only an ~8% contribution of tracks from the Forward 
Track Detector. The HCM current hemisphere, in contrast, has an energy scale of W/2 
rather than Q/2 and generally can not be seen in its entirety. This makes comparisons 
with a complete e+e~ interaction hemisphere somewhat easier in the Breit frame. 

When QCD corrections to the parton model are considered, the incident parton carries, 
in general, a proportion of the incident proton's momentum that is larger than x and has 
an energy larger than Q/2. Furthermore, since the four-momentum of the photon is fixed 
the energy seen in the current hemisphere may be either greater than or less than Q/2. 
In an e^e~ interaction it is also the case that the energy in any given hemisphere may 
be greater than or less than E* /2. The effects of final state radiation in ep DIS and e^e^ 
interactions are similar, but other LO QCD processes which affect DIS, initial state QCD 
radiation and boson gluon fusion (BGF), have no equivalent in e^e~ annihilation. Note 
that, to order a^, i.e. two exiting partons as well as the spectator system, DIS events 
may have no energy in the current hemisphere when the incident parton has energy above 
Q, other than that due to hadronisation P]. Momentum conservation ensures that this 
situation never exists in any given hemisphere of an e^e~ interaction. 



5 Evolution of the Fragmentation Function 

The ep Breit frame equivalent of the e+e" scaled hadron momentum Xp = 2phadron/E* is 
Xp = 2phadron/ Qi where only hadrons in the current hemisphere are considered. The event 
normalised distribution D{xp,Q'^) = (l/Nevts) x dn^^^^f^^/dxp, the fragmentation function, 
characterises the processes by which partons shower and then hadronise. In this paper, 
the intention is to present the spectra of charged particles originating from the primary 
vertex following the fragmentation of light quarks. 

At high enough energies and for light enough final state particles, fragmentation func- 
tions approximately scale and are "soft", rising rapidly at small Xp and peaking near 
Xp = 0. As the energy of the initial parton increases, D{xp, Q"^) evolves into an even softer 
function with increased probability for low Xp hadrons at the expense of high Xp. This 
scaling violation can be seen with the low and high Q^ data as a function of Xp in Fig. |I](a) 
or in Fig. |^ where the fragmentation function is plotted as a function of Q for different 
intervals of Xp. Typically, in this analysis, the intervals of Xp are much greater than the 
resolution (~6% rms). In these figures the HI data are compared with equivalent data 



from e~^e~ experiments plotted as a function of E*. Most published e~^e~ results 
refer to full event multiplicities. Here, as elsewhere in this analysis, track multiplicity 
data from e"'"e~ experiments have been halved to correspond to the fragmentation of 
one timelike quark. Note that HI as a single experiment is able to measure the violation 
for a wide range of Q. Fig. |^ also displays a prediction from a leading order QCD plus 
hadronisation Monte Carlo (DJANG06 [jT^) calculation. The scaling violation effect has 
a similar origin to the scaling violations in structure functions, and with next to leading 
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order calculations may be used as a test of perturbative QCD [|T^. 
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Figure 1: The current hemisphere fragmentation function as a function of (a) Xp and (b) 
^ shown separately for the low (open circles) and high (closed circles) data samples. 
Statistical and systematic errors are added in quadrature. The dashed and dotted lines 
show the results of Gaussian fits. 

When the fragmentation function is plotted as a function of the variable ^ = ln(l/a;p) 
the turn-over region is expanded. According to the expectations of the Modified Leading 
Logarithmic Approximation (MLLA) and Local Parton Hadron Duality (LPHD), D{C,) is 
Gaussian in the neighbourhood of the peak ||2^ . Even the high statistics of this, compared 



with our earlier, analysis see no significant deviation from this Gaussian behaviour, as is 
clear in Fig. |l](b). Most track measurement problems occur with particles of very high or 
very low momenta and these are correlated with low and high ^ respectively. To avoid 
any dependence on these areas the fits of this analysis are restricted to a region within 
one unit of the dimensionless ^ on either side of the mean, but tests in which this range 
is varied by ±20% show no systematic difference. 

The evolution of the fragmentation function may be summarised by the Q, or E*, 
dependence of the peak and width (dispersion) values of the fitted Gaussian as is shown 
in Fig. and in Table |I[ It is notable that the systematic effects discussed in section 3 
contribute only a total 2% uncertainty to the peak value and have no significant effect on 
the width measurement. The results are compatible with those published earlier, which 
had considerably lower precision, as well as with those of an analysis by the ZEUS 
collaboration. The results are also compatible with various e^e~ experiments pT|, |2^, 



where the relevant evolution variable is the centre of mass energy, E*. For this comparison 
and for the consideration of multiplicities in the next section, these published e~^e~ data 
have been re-fitted by us to be in a directly comparable form. As well as the previously 
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Figure 2: Scaling violations of the fragmentation function as a function of E* for e^e^ 
results (starred symbols), and as a function of Q for HI results (all other symbols) each 
within the same indicated intervals of the scaled momentum, Xp. Note that the data for 
the three lowest intervals of Xp are multiplied by factors of 10, 5 and 3 respectively for 
clarity. Statistical and systematic errors are added in quadrature. The solid lines show 
the prediction of a DJANGO Monte Carlo calculation. 



mentioned factor of two to convert e^e^ event multiplicities to quark multiplicities, we 
also reduce published total multiplicity results by 8% to account for K'^ and A decay 
tracks |2^, 24 1, and by a variable factor below 3% to correct for the increased multiplicity 
from h quark fragmentation ||25| in e^e~ . 

Assuming gluon coherence MLLA/LPHD predicts a dependence of the evolution of 
the peak and width of the distribution on the dimensionless variable Y = \n{Q / A ^jj) , 
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* 2 




Q,E (GeV) 

Figure 3: HI results (solid symbols, with statistical and systematic errors combined in 
quadrature) showing the evolution of (a) the peak and (b) the width of the fragmentation 
function as a function of Q compared with e~^e~ results (open symbols) as a function of 
the centre of mass energy, E* . The solid line is a fit to MLLA/LPHD expectations. 



where A ^jj is a scale parameter determining a presumed cut-off of a parton shower. The 
prediction ||20| gives the peak, ^peak and the width, ^width, of the Gaussian approximation 
to be 

^peak = 0.5Y + C2VY + IC 



e 



width 



Y-2/2CI 



where ci = ^36Njb and Ca = B^^b/N^ with b = ^N^ - |A^/ and B -- 
^Nf /N^)/b are constants dependent only on the number of contributing colours, Nc, and 
flavours, Nf within the parton shower. For ease of comparison, we follow |22| in the 
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(Q) 

(GeV) 


Total Current Hemisphere 


Energy Flow Selected 


{n) 


reak 


Width 


(n) 


reak 


Width 


3.63 


1.18±0.07 


1.52±0.05 


0.80±0.06 


1.77±0.11 


1.51±0.05 


0.76±0.07 


4.12 


1.38±0.08 


1.53±0.04 


0.74±0.05 


1.85±0.11 


1.51±0.03 


0.68±0.04 


5.23 


1.62±0.08 


1.72±0.03 


0.75±0.03 


O "1 O 1 1 1 

2.13±0.11 


1.70±0.03 


0.76±0.03 


6.90 


2.19±0.12 


1.84±0.04 


0.78±0.04 


2.60±0.15 


1.83±0.04 


0.77±0.05 


8.20 


2.59±0.16 


1.97±0.08 


0.90±0.10 


3.10±0.20 


1.96±0.07 


0.87±0.10 


9.32 


3.16±0.25 


2.04±0.10 


0.87±0.12 


3.49±0.29 


2.06±0.09 


0.84±0.11 


12.5 


3.65±0.32 


2.26±0.11 


1.09±0.18 


3.93±0.35 


2.23±0.10 


1.03±0.17 


14.6 


4.13±0.34 


2.43±0.08 


1.05±0.13 


4.45±0.30 


2.41±0.07 


1.04±0.13 


17.9 


4.05±0.34 


2.51±0.06 


0.97±0.09 


4.28±0.36 


2.49±0.06 


0.98±0.10 


25.0 


5.23±0.43 


2.77±0.08 


1.08±0.14 


5.28±0.44 


2.80±0.10 


1.18±0.18 


41.2 


7.32±0.65 


3.02±0.14 


1.26±0.23 


7.19±0.65 


2.95±0.14 


1.19±0.20 



Table 1: The average charged multiphcity, and the peak (^peak) and width (^width) of the 
fragmentation function as a function of Q^, for the total current hemisphere of the Breit 
frame both with and without the energy flow selection discussed in section 6. The errors 
are the sum of statistical and systematic uncertainties in quadrature. 



assumption that this will be dominated by the light quarks and set Nj = 3. The term 
/C contains higher order corrections and is expected to be roughly constant and of order 
1. The result of a simultaneous fit (x^/NDF of 10/20) to the peak and width values 
obtained from the present HI data alone is shown as a solid line in Fig. ^ and yields 
the values A^jj^ = 0.21 ± 0.02 GeV and JC = —0.43 ± 0.06, in agreement with the values 
A^jj = 0.21 ± 0.02 GeV and /C = —0.32 ± 0.06 obtained in an analysis of the ^peak 
evolution of combined e"'"e~ data. 



6 The Average Charged Particle Multiphcity 

The area underneath the fragmentation function is the average charged multiplicity and 
this is given, after all corrections, in Table |I|. It is shown as a function of Q in Fig. ^, with 
ZEUS results and with a curve which accurately parameterises the average charged 
multiplicity from many e~^e~ experiments as a function of E*. The error associated with 
this fit at any given energy is at the percent level. 

There are predictions that both the absolute average charged DIS multiplicity 
and the shape of the evolution with should agree with that of e~^e~ experiments. 
As previous results indicated, although there was very good agreement between HI and 
ZEUS [|] analyses and those of e~^e~ at high Q{E*), the DIS analyses give smaller average 
charged multiplicities at low values of Q{E*). The new data confirm that this is a large 
effect. 

The discussion of section 4 indicated that the LO processes which are present in ep 
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1 10 . 10^ 

Q,E GeV 

Figure 4: Average charged multiplicity in the current hemisphere of the Breit frame as 
a function of Q for this analysis (solid circles) and for the ZEUS 1993 analysis (open 
squares). Statistical errors and systematic errors are added in quadrature. The curve is 
a fit to many e~^e~ results as a function of the centre of mass energy, E* . 

but not in e~^e~ interactions can produce a depletion of the current region of the Breit 
frame. Pure QPM interactions or events with only final state radiation (i.e. e~^e~ like 
events) might both be expected to have a current region hadronic state characteristic of 
a quark recoiling with energy ~ (5/2 and with only a z component of momentum. BGF 
events with a depleted or empty current hemisphere or events with significant initial state 
QCD radiation would, on the contrary, have a current region hadronic final state with 
significantly less energy at an angle to the z axis. However, hadronisation effects might 
well mask these naive expectations. 

In order to investigate this experimentally, we add the four-momenta of all calorimeter 
energy clusters in the current region of the Breit frame. The energy component of the 
resultant four-momentum vector is defined to be Ez<q and the angle it makes with the z 
axis to be Qbf- We then plot Ez<o/Q against cosG^f, as shown in Fig. ^(a,b) for the 
low and high samples. From the above discussion we would expect that e^e~ like 
events might cluster near (cosB^i;' = —l.E^^^o/Q = 1/2) in this plot, whereas events 
from the other LO processes would either not appear in this plot or be spread away from 
this point, typically towards lower energies and larger angles in the current hemisphere, 
that is to the lower right-hand corner of this plot. The high sample, where as is small. 
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X is large, and the proton has a much smaller gluon content, shows obvious clustering 
near the (-1,0.5) point, while this is much less evident at low . 

Monte Carlo simulations qualitatively confirm the expectations but it is hard to be 
quantitative, both because of the difficulty in strict definitions of LO QCD event classes 
and because experimental resolution has a large contribution to event migration on this 
plot. Nevertheless, it is possible to investigate the effect of choosing ever more e^e~ 
like events by selecting according to the variable R, where = {Ez<:o/Q — 0.5)^ + 
{cosQbf + 1)^- In order to compare with e~^e~ results, we divide {n{R)), the average 
multiplicity of events in an annulus of width AR = 0.1 at radius R, by the equivalent 
single-hemisphere e~^e~ multiplicity. There is a very small variation of (Q) with R, so 
this e~^e~ multiplicity is evaluated at an energy corresponding to the average value of Q 
of the HI data in that annulus. The behaviour of this ratio as a function of R can be 
seen, with statistical error bars only, in Fig. |^(a,b). The systematic errors corresponding 
to the shaded region are as discussed in section 3. Events with a totally empty current 
hemisphere have zero multiplicity but no defined value of R, and are thus not included in 
these plots. As i? — >■ the average charged hadronic multiplicity increases, possibly even 
beyond the corresponding e~^e~ values, although it is not inconsistent with e~^e~ values 
given the systematic errors involved. Our Monte Carlo simulations also show this general 
behaviour. 
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Figure 5: The total energy oj the summed calorimeter cluster four-momentum vectors in 
the current hemisphere of the Breit frame is plotted as a fraction of the event Q against 
the polar angle of the resultant vector, for (a) the low and (b) the high data sample. 
The annuli define the R variable for Fig. |^ and the solid line indicates the Breit frame 
energy selection referred to in the text. 



The above result and discussions motivate a cut to remove extreme low energy and 
off-axis events. This is shown as a solid line in Fig. |^(a,b) joining the points (-1,0) and 
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Figure 6: The HI average charged multiplicity as a function of the radius defined in the 
text normalised to that of the single hemisphere expectation for e+e" events for (a) the 
low and (h) the high data sample. The HI points have statistical error bars only. 
The shaded regions show the total systematic error. 



(0,0.5) in this plane and is termed the 'Breit frame energy flow selection' in the rest 
of this paper. It produces samples of 11 735 events and 1 005 events at low and high 
respectively, which are less affected by the low-order QCD processes peculiar to DIS. This 
selection is by no means as exact as the technique of letting — > but in circumstances 
where the data are already binned in other variables, e.g. energy spectra or multiplicity 
distributions, gives an indication of whether or not distributions move in the direction of 
closer similarity between this data and e"'"e~ data. 

As is clear from Table |1|, the peak and width values of the fragmentation function are 
insensitive to this energy flow selection procedure. The average multiplicity, however, is 
sensitive and as expected increases with this selection, especially at low Q"^. 



7 Invariant Energy Spectra 

The evolution of the fragmentation function at low energies, or high values of ^, may be 
better studied using the event-normalised "invariant" spectrum, {1/ Nevts)^Edn%^^j^g/ d^p. 
According to the predictions of perturbative QCD based on the MLLA and LPHD, the 
hadronic spectrum at low momentum should be nearly independent of the energy of 



the parent parton Here we follow in taking d^p = ^T^Phadron^Vhadron and by 
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computing the energy, E, of each track as if it had a mass Qo{=270 MeV~ Aqcd) given 



by a presumed cut-off mass for the parton shower, i.e. E = ^Ql + Vhadron- The data are 
corrected for acceptance with the same bin- by-bin method described in section 3. The 
correction function is both smooth and close to unity, having an average value of around 
1.15. Although there is always a loss of track acceptance at low laboratory momentum 
the boost reduces this loss in the Breit frame. This means that ep experiments can get 
closer to the critical low energy limit without serious loss of acceptance, than can e^e~ 
experiments. 

The results are shown in Fig. 0(a,b) where the indicated errors contain contributions 
from statistical and (negligible) momentum measurement effects as well as systematic 
effects that have already been discussed and which grow slightly with track energy. The 
solid lines show MLLA/LPHD predictions at the relevant Q values for a simple quark 
behaviour using the calculations of I^H] which the authors show are in agreement with 
e^e~ data. It is clear that at low these predictions for the behaviour of quark fragments 
are in disagreement with the data but are significantly closer to data subject to the Breit 
frame energy flow selection described in the previous section. This selection makes little 
difference at high Q^, where agreement is better anyway. The QCD predictions are very 
sensitive to the 'running' of the strong coupling constant, a^, as the energy scale alters 
in the parton shower. A prediction is shown for the same average energy, normalised to 
interpolated {n^) multiplicities in e^e~ experiments, but utilising a fixed value of as (set 
to 0.215, a value which roughly describes the energy dependence of (jt^) and the slope 
of the first moment of the energy spectrum in e^e~ experiments [^). The data show 
that such an assumption is clearly untenable. This calculation is only in leading order so 
this cannot yet be turned around to compute a meaningful value of Aqcd- Data subject 
to the energy flow selection are re-plotted in six intervals of in Fig. |^(c) in order to 
display the clear evidence for a common limit at low particle energy independent of event 
momentum transfer. 



8 Charged Particle Multiplicity Distributions 

Charged hadron multiplicity distributions are considerably more difficult to measure than 
average values. In this analysis we have tested a number of ways of extracting true 
distributions from observed distributions. The best method is found to be a simple bin-by- 
bin correction procedure on the multiplicity distribution. This has, however, the problem 
of not modelling correlations which then gives the method a dependence on the Monte 
Carlo generator. 

We tried two further methods to solve this difficulty. It is straightforward to use 
a Monte Carlo simulation to form a generated, or 'true', to reconstructed distribution 
transfer matrix, which is physics- generator independent. Inversion of this matrix in order 
to extract a true distribution from a raw input distribution has many problems due to 
the small number of high multiplicity events, resulting in the matrix becoming singular 



and unstable. We have also tried a method [T^, pO[ of approximating the inverted matrix 



with its inverse and then re-weighting the input Monte Carlo distribution in successive 
approximation iterations to minimise dependence on the initial generated distribution. 
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Figure 7: The invariant charged hadron energy spectrum in the current hemisphere (a) 
at low and (b) at high Q^. The data for all events are shown as open circles and 
those utilising the Breit frame energy flow selection as solid circles. The solid line is 
the prediction of MLLA/LPHD and the dashed line is the corresponding expectation for a 
non-running coupling constant, ag. In (c) the low data are subjected to the energy flow 
selection hut subdivided into the six indicated intervals of and plotted with a regular 
incremental spacing of 0.5 GeV on the abscissa. The solid line MLLA/LPHD expectation 
is calculated at the (Q) for each distribution. The error bars show the sum of statistical 
and systematic errors added in quadrature. 
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The systematic errors of all three techniques have been estimated by degrading the 
output distributions back to the 'raw' (reconstructed) data level using the rehable 'gen- 
erated to reconstructed' matrix. This can then be compared with the original, real input 
distribution. We find that the larger dependence on the generated distribution of the bin- 
by-bin correction method is still offset by overall smaller systematic errors and the ability 
to get meaningful corrections at high multiplicity when the statistics are low. We thus 
choose this method to present our results. Quantitatively, the average correction made 
per bin is 16%(32%) at low (high) with an error of method assessed as 10%(17%). 

The resulting multiplicity distributions are given in Table |^ and shown in Fig. ^ in the 
form of P{n), the probability of observing a charged multiplicity n. These results are for 
the multiplicity distribution in the current hemisphere of all events and do not include 
any energy flow selection. The bins of x and are exclusive but, because of the event 
selections, are not totally occupied at the edges. This leads to only minor discrepancies in 
the average values of x and between neighbouring bins. In general the observation is 
that the results are compatible with a gradual shift to higher multiplicities with increasing 

but with only a small change as a function of x which might be due to minor bin edge 
effects. It is interesting to note the many zero-charged multiplicity events, especially at 
low where we might expect higher-order QCD processes to dominate. 

Integrating over x gives smaller statistical errors and a greater reach in multiplic- 
ity. The evolution with Q of the shape of the multiplicity distribution may be studied 
by normalising with the average charged hadron multiplicity to give the KNO form 
"^{z) = (n) X P{n), where z = n/ (n). Equivalent single-hemisphere e~^e~ distrib- 
utions [0, 0, ^ exhibit clear scaling in this variable but only for high {E* > 20 GeV) 
energies. In an analysis [|T7| of multiplicity distributions in a limited region of the hadronic 
centre of mass frame, the HI collaboration has also observed good agreement with a com- 
mon distribution as W increases between 80 and 220 GeV. The results of this analysis in 
the Breit frame, shown in Fig. §(a), are for Q values much lower than the W values and 
show a clear violation of KNO scaling, as also reported by the ZEUS collaboration 
For reference, we display two sets of lines joining e~^e~ data at an energy close to our 
largest Q range and also at LEP energies ||2^. As with e~^e~ data, the HI low 
distribution starts with a flatter curve and a large contribution from low multiplicities. 
Then, as increases, the data get closer to the high energy scaling e~^e~ distributions. 

In Fig. |^(b) the distribution for low events is compared with results from the e~^e~ 
JETSET parton shower model, using parameter settings determined by the DELPHI 
collaboration at an energy set to be close to the mean value of the low data. The 
data distribution is again much flatter but a better agreement is achieved when using the 
Breit frame energy flow selection discussed in section 6. This KNO scaling violation is 
also seen in ep DIS Monte Carlo generated events, where it is to a large extent caused 
by leading order perturbative QCD events, especially BGF events, which, as has already 
been noted, have signiflcantly lower multiplicities in the current hemisphere of the Breit 
frame. It thus appears that ep and e~^e~ data share a flattening of the KNO distribution 
at low energies but that there is an additional sensitivity to the effects of higher order 
QCD processes in the Breit frame which do not contribute to e~^e~ events. 
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Figure 8: The event-normalised charged hadron multiplicity distribution in the current re- 
gion of the Breit frame, with statistical and systematic errors added in quadrature. Mov- 
ing from bottom to top the distributions are for 12 < < 30 GeV^, 30 < < 80 
GeV^ and 100 < < 500 GeV^, and from left to right the ranges in Bjorken-x are 
6 ■ 10"^ < X < 2 ■ 10-'\ 2 ■ IQ-^ <x <1- 10-2 and 1 ■ 10"^ < x < 2 ■ 10"^ The average 
value of Q and x for each interval is given in Table |^. The average multiplicity and dis- 
persion are shown at the top of each distribution. The histograms show the prediction of 
a D J AN GO Monte Garlo calculation. 
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(Q) 

(GeV) 


<x> 
(xlO^) 


P(0) 


P(l) 


P(2) 


P(3) 


P(4) 


P(5) 


P(6) 


P(7) 


P(8) 


P(9) 


P(10) 


P(ll) 


4.15 


1.26 


0.312 ±0.025 


0.309 ±0.016 


0.219 ±0.027 


0.116 ±0.009 


0.031 ± 0.005 


0.005 ± 0.001 


0.0012 ± 0.0006 


- 


- 


- 


- 


- 


4.53 


2.55 


0.246 ± 0.070 


0.281 ± 0.025 


0.233 ± 0.042 


0.133 ±0.031 


0.055 ±0.011 


0.030 ±0.012 


0.0012 ± 0.0020 












6.20 


1.58 


0.194 ±0.023 


0.250 ±0.026 


0.205 ±0.025 


0.182 ±0.023 


0.079 ±0.017 


0.036 ±0.011 


0.031 ±0.018 












7.04 


4.10 


0.150 ±0.031 


0.195 ±0.013 


0.229 ±0.025 


0.199 ±0.020 


0.116 ±0.026 


0.049 ±0.012 


0.023 ± 0.006 


0.007 ±0.004 


0.003 ±0.003 








14.5 


6.82 


0.041 ±0.009 


0.083 ±0.019 


0.119 ±0.024 


0.236 ± 0.049 


0.169 ±0.024 


0.124 ±0.024 


0.106 ±0.027 


0.035 ± 0.017 


0.030 ±0.017 


0.017 ±0.011 






17.3 


2.34 


0.020 ±0.011 


0.052 ±0.014 


0.093 ±0.018 


0.171 ±0.031 


0.184 ±0.026 


0.151 ±0.032 


0.135 ±0.029 


0.075 ±0.026 


0.048 ± 0.022 


0.021 ±0.010 


0.018 ±0.013 


0.007 ±0.005 



Table 2: The probability per event, P{n), of observing a charged hadronic multiphcity 
of n in the current hemisphere of the Breit frame as a function of Q and x. The errors 
are the sum of statistical and systematic errors in quadrature. 
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Figure 9: (a) The charged hadronic multiplicity distribution for HI data expressed in terms 
of the KNO variables for 12 < < 30 GeV^ (open circles), 30 < < gO GeV^ (squares) 
and 100 < < 500 GeV^ (solid circles), (b) A repeat of the 12 < < 30 GeV^ HI 
data with (solid triangles) and without (open circles) the Breit frame energy flow selection. 
Statistical and point to point systematic error bars are added in quadrature in both figures. 
The dotted and solid lines represent the single hemisphere KNO functions for e'^e~ data 
at E* —I4 and 91 GeV respectively and the dashed line shows the result from the e'^e" 
JETSET Monte Garlo at an E* value close to the Q value of the ep data. 
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9 Summary 



The Breit frame fragmentation function analysis has been repeated with the higher statis- 
tics of 1994 data and is now extended to give multiphcity distributions as functions of 
both X and Q^. Invariant energy spectra are shown to be incompatible with a calculation 
utilising a constant value for as- We discuss a method of discriminating against events in 
kinematic regions dominated by QCD processes which do not contribute to e'^e" interac- 
tions. We show that some discrepancies between ep Breit frame current hemisphere and 
e~^e~ average charged multiplicities, including an observed violation of KNO scaling at 
low Q^, decrease as conditions are progressively made more comparable. With a selection 
which gives similar kinematic conditions to that of a single hemisphere e~^e~ annihilation, 
the ep energy spectra can be shown to be compatible with a common, Q^-independent, 
low momentum limit. In summary, the data give strong support to the concept of quark 
fragmentation universality. 
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